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Abstract 

This paper focuses on validating the proposed Neutrosophic 
Set (NS) approach of Magnetic Resonance Image (MRI) 
denoising based on structural similarity such as Structural 
Similarity Index (SSIM) and Quality Index based on Local 
Variance (QILV). The Neutrosophic Set approach of median 
filter is used to reduce the Rician noise in MR image. The 
experiments have conducted on real MR image with Rician 
noise added. The visual and the diagnostic quality of the 
denoised image is well preserved. The performance of this 
filter is compared with median filter and non local mean filter 
(NLM). 

1 Introduction 

Magnetic resonance imaging (MRI) is the most powerful 
imaging technique [22] developed to study the structural 
features and the functional characteristics of the internal body 
parts. The diagnostic and visual quality of the MR images are 
affected by the noise added while acquisition. This is 
problematic for further tasks such as segmentation of 
important features; classification of images for computer 
aided diagnostics, three dimensional image reconstruction 
and image registration. Therefore, denoising should be 
performed to improve the image quality for more accurate 
diagnosis. 

The MR image is commonly reconstructed by computing the 
inverse Discrete Fourier Transform of the raw data. The 
signal component of the measurement is present in both real 
and imaginary channels. Each of the orthogonal channels is 
affected by additive white Gaussian noise. The noise in the 
reconstructed complex valued data is thus complex white 
Gaussian noise. Most commonly, the magnitude of the 
reconstructed MRI image is used for visual inspection and for 
automatic computer analysis. Since the magnitude of the MRI 
signal is the square root of the sum of the squares of two 
independent Gaussian variables, it follows Rician distribution 
[8]. The Rician noise is signal dependent and is therefore 
difficult to separate from the signal. 

Numerous approaches of denoising MR images have been 



proposed including approaches based on anisotropic diffusion 
[7,10,19,20,26], the wavelet transform [2,6,15,16,17,25], 
bilateral and trilateral filtering [9,23,24], the non-local means 
algorithm [3-5,11-13] in the literatures. This paper is an 
extended work of [14]. The Neutrosophic Set approach of 
MRI denoising is validated in terms of structural similarity 
(SSIM) [21] index and quality index based on local variance 
(QILV) [1]. As already mentioned in [1, 21], peak signal to 
noise ratio (PSNR) is not the reliable quality metric with 
respect to human visual system (HVS). SSIM is an objective 
quality measure based on the structural content of the image 
that engulfs the properties of HVS. Medical images contain 
delicate structural features that contribute more to the 
diagnosis. QILV will provide the statistical similarity between 
the denoised and the original image. 

2 Neutrosophic MRI denoising 

Neutrosophy, a branch of philosophy introduced in [18] as a 
generalization of dialectics, studies the origin, nature and 
scope of neutralities, as well as their interactions with 
different ideational spectra. Neutrosophy theory considers 
proposition, theory, event, concept or entity, <A> is in 
relation to its opposite <Anti-A> and the <Neut-A> which is 
neither <A> nor <Anti-A>. The neutrosophy is the basis of 
the neutrosophic logic, neutrosophic probability, and 
neutrosophic set and neutrosophic statistics [18]. In 
neutrosophic set, the indeterminacy is quantified explicitly 
and the truth-membership, indeterminacy-membership and 
falsity-membership are independent. The neutrosophic set is a 
general formal frame work which generalizes the concept of 
the classic set, fuzzy set, interval valued fuzzy set, 
intuitionistic fuzzy set, and interval valued intuitionistic fuzzy 
set, paraconsistent set, dialetheist set, paradoxist set and 
tautological set [18]. The neutrosophic set had been applied 
into image processing such as denoising [27], segmentation 
[28, 29]. The definition of a neutrosophic set and its 
properties are described briefly. 

2.1 Neutrophic Set 

Definition 1 (Neutrosophic Set). Let U be a Universe of 
discourse and a neutrosophic set A is included in U . An 
element x in set A is noted as x(T,I,F). T.fFdXQ real 
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standard and non standard sets of ]-(),l + [with Definition 3 (Neutrosophic image entropy). Neutrosophic 

entropy of an image is defined as the summation of the 
sup T = t_ sup, inf T = t _ inf, sup I = i _ sup, inf I = i _ mf, entropies of three subsets T IandF: 



supF /_ sup, inf F f inf and z?_sup ^_sup + 
i _ sup + / _ sup, n _ inf = t _ inf + i _ inf + / _ inf . 

T, I and F are called the neutrosophic components. The 
element x(T,I,F) belongs to A in the following way. It is 
t% true in the set, z’% indeterminate in the set, and /% false 
in the set, where t varies in T , i varies in / and / varies 
in F . 

2.2 Transform the image into neutrosophic set 

Definition 2 (Neutrosophic image). Let U be a Universe of 
discourse and FT is a set of U , which is composed by bright 
pixels. A neutrosophic image P NS is characterized by three 
membership sets T 7 , /, F . a pixel P in the image is described 

as F(r,/,F) and belongs to W in the following way: It is 
t true in the set, z indeterminate in the set, and f false in the 
set, where t varies in F , i varies in / and / varies in F . 
Then the pixel F(z , /) in the image domain is transformed 
into the neutrosophic set 

domain P NS (/, j) = {T (/, /), I (/, j), F(i, j)}.T (/', j ), /(/, /) 
and F(z,y) are the probabilities belong to white pixels set, 
indeterminate set and non white pixels set respectively, which 
are defined as: 



T(i,j)= g ( I,J ^ gmin 
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S(i, j ) = abs(g(i, j ) - g(i, j)) 


(4) 


F(i,j) = \-T(i,j) 


(5) 



where g(z, j) is the local mean value of the pixels of the 
window. S(i,j) is the absolute value of difference between 
intensity g(z, /) and its local mean value g(i, j) . 

2.3 Neutrosophic image entropy 

For a gray image, the entropy is utilized to evaluate the 
distribution of the gray levels. If the entropy is the maximum, 
the intensities have equal probability. If the entropy is small, 
the intensity distribution is non-uniform. 



En NS = En T + Fz? 7 + En F 
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where En T , Enj and En F are the entropies of sets 
T 7 , / and F respectively. p T (i) , pj (i) and p F (i) are the 
probabilities of elements in T 7 ,/ and F respectively, whose 
values equal to z . 

2.4 y - median filtering operation 

The values of /(z, j) is employed to measure the 
indeterminate degree of element F^ (/,/) . To make the set 

/correlated with T and F, the changes in T and F 
influence the distribution of element in / and vary the 
entropy of / . 

Definition 4 ( y - median filtering operation). A y — 

median filtering operation for P NS ^Ns(r) , is defined as: 



Adr)=7>(f( r ), />),%)) 


(10) 
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F y (z, y) = median { F(m , /?)} 
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where 8^{i,j) is the absolute value of difference between 

intensity T (/, /) and its local mean value f (/, /) at 
(/, /) after ^ - median filtering operation. 

The summary of neutrosophic set approach of MRI denoising 
is described as below (see figure 1): 



Step 1: Transform the image into NS domain; 

Step 2: Use y - median filtering operation on the true subset 

T to obtain f y ; 

Step 3: Compute the entropy of the indeterminate subset 
I y ,En- (/) ; 

En j (i + 1) - En- (i) 

Step 4: if < 8 , go to Step 5; 

En } (0 

1 r 

Else T = T y , go to Step 2; 

Step 5: Transform subset f y from the neutrosophic domain 
into the gray level domain. 



Magnetic Resonance Image 



Add Uniform Rician Noise 





Figure 1: Neutrosophic Set approach of MRI denoising 



3 Results and Validation 

Two brain MR images of 5 th slice, T2 Axial with 
repetition/ echo time (TR/TE) 9000/87 ms and T1 Sagittal 
with TR/TE 552/17 ms which are obtained using SIEMENS 
1.5 T MRI scanner for a 2 years old male child are considered 
for the experiments. This method is compared with median 
filter and Non Local Mean (NLM) filter proposed by Buades 
et al [3]. The results of the different denoising filters on a 
brain MRI corrupted with Rician noise is shown in Figure 2. 




(e) 



Figure 2: Results of different denoising filters on a brain 
MRI corrupted with Rician noise of standard deviation a n =15. 
T2 Axial (left) and T1 Sagittal (right) slices. From top to 
bottom: (a) the original image, (b) the noise corrupted image, 
and the results of (c) median filter, (d) the non local mean 
filter and (e) the neutrosophic set filter. 
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The performance of the denoising algorithm is measured by 
using the peak-signal-to-noise ratio (PSNR), the Structural 
Similarity (SSIM) index] and the Quality Index based on 
Local Variance (QILV). The PSNR will give the pixel by 
pixel similarity. The peak signal to noise ratio in decibel (dB) 
is measured using the following formula: 



PSNR = - 10 log 






HxWx 255 z 



(18) 



where /(z, j) and I d (z,y) represent the intensities of pixels 
(z, j ) in the original image and denoised image respectively. 
The higher the PSNR, the better the denoising algorithm is. 



SSIM and QILV give the measure of the structural similarity 
between the original and the denoised images and are in the 
range of 0 to 1. The SSIM works as follows: Let x and y be 



two non negative images, where as one has perfect quality. 
Then, the SSIM can serve as a quantitative measure of the 
similarity of the second image. The system separates the task 
of similarity measurement into three comparisons: luminance, 
contrast and structure. It can be defined as 
, x ( 2jU xJ u + C x )(2cr + C 2 ) 

SSIM(x,y) = ( 0 V V _ \(\ o A (19) 



\Mx + Ml + C 



ft 



<jI +<j 2 y + c 2 



where ju x and fu y are the estimated mean intensity and <j x and 



G y are the standard deviations respectively. a xy can be 
estimated as 

l N , v 

°xy = T7— 7 Z ( x i - Mx ) L " My j (20) 

7V 1 7=1 

C 1 and C 2 in eqn. 19 are constants and the values are given as 



C, ={K ] L) 1 and C 2 =(k 2 L f where K U K 2 « 1 is a small 
constant and L is the dynamic range of the pixel values (255 
for 8 bit gray scale images). 

The QILV between two images x and y can be defined as 



QILV(x,y) 



^MvxMVy 

2 2 
MVx + MVy 



^ G Vx G Vy G VxVy 

G Vx + G Vy <J Vx <J Vy 



( 21 ) 



where /u Vx and n Vy are the estimated means of the local 



variance and g Vx and a Vy are the standard deviations of the 



local variance respectively. G VxVy is the covariance between 

the two images. The performance comparison of Median, 
NLM, and NS Median Filters based on the PSNR in dB, 
SSIM and QILV values for T2 Axial and T1 Sagittal slices of 
brain MRI corrupted with different levels of Rician noise are 
given in Figure 3 and Figure 4 respectively. The numerical 
values of the quality measures: PSNR in dB, SSIM and QILV 
are listed in Table 1 for the different denoising filters. As 
already discussed in [14], According to the values of PSNR, 
at High SNR (low noise level), NLM performs well than NS 
filter and at low SNR (high noise level), NS filter performs 
better than that of NLM. This result is justified by the quality 
measures SSIM and QILV. 







(d) 



Figure 3: Comparison of Median, NLM, and NS Median 

Filters based on the PSNR in dB values: T2 Axial (a) High 
SNR (b) Low SNR; T1 Sagittal (c) High SNR (d) Low SNR. 



4 






Image 


Noise 


PSNR 


SSIM 


QILV 






Median 


NLM 


NS Median 


Median 


NLM 


NS Median 


Median 


NLM 


NS Median 




2 


23.18 


35.65 


32.69 


1 


1 


1 


0.9721 


0.9947 


0.9940 


T2 

Axial 


10 


23.14 


34.36 


30.41 


0.8920 


0.9805 


0.9564 


0.9048 


0.9921 


0.9784 


20 


22.79 


28.36 


26.84 


0.8499 


0.9328 


0.9096 


0.8074 


0.9716 


0.9644 


30 


22.04 


24.81 


24.19 


0.7911 


0.8644 


0.8481 


0.7040 


0.9309 


0.9226 




40 


21.09 


21.30 


21.99 


0.7199 


0.7819 


0.8145 


0.6000 


0.8685 


0.8909 




50 


19.94 


20.27 


21.05 


0.6459 


0.6949 


0.7522 


0.5109 


0.7749 


0.8133 




60 


18.81 


19.64 


20.02 


0.5732 


0.6092 


0.6961 


0.4431 


0.6766 


0.7612 




70 


17.74 


18.25 


19.90 


0.5081 


0.5336 


0.6352 


0.3913 


0.5730 


0.6879 




2 


27.75 


46.56 


33.40 


1 


1 


1 


0.8562 


0.9985 


0.9970 




10 


27.24 


35.23 


31.42 


0.7537 


0.9756 


0.8963 


0.7449 


0.9883 


0.9665 


T1 


20 


25.89 


29.19 


28.14 


0.7037 


0.9223 


0.8429 


0.6384 


0.9694 


0.9469 


Sagittal 


30 


24.20 


25.58 


25.59 


0.6440 


0.8457 


0.8244 


0.5500 


0.9317 


0.8909 




40 


22.48 


22.90 


23.28 


0.5769 


0.7586 


0.7670 


0.4814 


0.8734 


0.8593 




50 


20.86 


20.80 


22.00 


0.5149 


0.6700 


0.6998 


0.4365 


0.8038 


0.8313 




60 


19.44 


19.07 


20.36 


0.4517 


0.5842 


0.6359 


0.3984 


0.7615 


0.7643 




70 


18.08 


17.55 


19.11 


0.3985 


0.5078 


0.5795 


0.3826 


0.6489 


0.7370 



Table 1 : Performance Comparison of Median, NLM and NS Median Filters 





(a) (b) 





Figure 4: Comparison of Median, NLM, and NS Median Filters based on the quality measures SSIM and QILV: T2 Axial (a) 
SSIM (b) QILV; T1 Sagittal (c) SSIM (d) QILV. 



5 



That is, According to the structural similarity measures, by 
applying the Neutrosophic Set approach of median filter, the 
detailed structures are well preserved for low SNR compared 
with Median filter and Non local mean filter. 

4 Conclusion 

The proposed denoising technique based on Neutrosophic Set 
for reducing Rician noise from MR image have been 
validated in this article by using the structural similarity 
measures (SSIM, QILV). Experiments have been carried out 
on real MR images. This filtering method tends to produce 
good denoised image not only in terms of visual perception 
but also in terms of the quality measures such as PSNR, SSIM 
and QILV. This filter performs better than Median filtering 
method for reducing the Rician noise with different noise 
levels. Further, and also it outperforms the Non Local Mean 
approach when the noise level is high (low SNR). 
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